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Abstract 

We study the motion of elastic networks driven over a random substrate. Our 
model which includes local friction forces leads to complex dynamical behavior. We 
find a transition to a sliding state which belongs to a new universality class. The 
phase diagram comprises of a pinned state, a stick-slip motion phase, and a free 
motion phase. 



1 Introduction 



The dynamics of elastic objects driven through an environment with random 
pinning interactions is complex and its present understanding still incomplete. 
Processes of this type occur in microscopic systems such as charge density 
waves driven by an electric field [1,2], superconductors subject to external 
magnetic fields [3], and polymers moving through a inhomogeneous environ- 
ment. Related processes on a macroscopic scale occur in the motion of geolog- 
ical faults [4-6], or in sliding friction between an elastic object and a rigid one 
[7]. Theoretical studies of these phenomena are usually modeled by an elastic 
network moving in a rigid environment or over a rigid substrate with random 
interactions between them. In the simplified case where one of the two objects 
is homogeneous, one may distinguish between random network models, where 
the random interactions are assigned to sites on the elastic network while 
the medium or substrate are homogeneous, and the random substrate, models 
where the random interactions are assigned to the sites of the substrate or the 
embedding medium while the elastic network is homogeneous. These two type 
of models are illustrated in Fig. 1 

Random network models were invoked to describe randomly pinned charge 
density waves (CDW) in solids [8]- [13] and were widely studied. Most recent 
work on random substrate models concentrated on the high velocity limit, 
where the distortions of the network induced by the random interactions are 
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mild [14]-[17], while the dynamical properties near the critical external field for 
depinning received less attention. We focus on the simplest discrete realization 

of the random substrate model containing no ad-hoc features such as veloc- 
ity dependent forces [5], and compare it to the analogous random network 
model. The study shows that the dynamics of the random substrate model 
above depinning threshold is interesting and complex, featuring several dis- 
tinct regimes. In contrast, in the random network case the elastic object follows 
a simple steady state motion. The pinned state of the random substrate model 
below threshold is characterized by a wide distribution of strain avalanches. 
Similar avalanches occur in the random network model, however the critical 
exponents describing the distributions of avalanche sizes are definitely differ- 
ent. 

The properties of the discrete random substrate model depend on three param- 
eters: The external field, the ratio of pinning strength to the elastic stiffness, 
and the characteristic size of the local pinning regions. We find that whereas 
in some parts of the relevant phase diagram the two models are approximately 
equivalent, through most of it there are pronounced differences between them, 
and the critical behavior belongs to distinct universality classes. 



2 The Model 

Sliding elastic objects are represented here by a discrete chain of particles 
connected by springs, (see Fig. 1) In the random network limit each particle 
is attributed a random pinning force, while the substrate is homogeneous. For 
an external force E switched on at t = acting uniformly on all particles, 
the dynamics of the displacements of the positions of the particles from their 
initial equilibrium positions: = Xi{t) — Xi{0) is described by the following 
overdamped equations of motion: 

random network 

dii/dt = max(0, + Ci-i - 2^0 + E- V{i)) (1) 

This realization of the random network model, also known as the ratchet model 
[9] successfully describes the dynamical phase transition of charge density 
waves from a pinned state to a DC current conducting state at a critical 
external field [9-12,18]. 

In the random substrate model the particles interact with the substrate via 
local pinning forces V{xi) , whereby the i'th particle at position Xi (in the 
substrate frame of reference) moves if the total force acting on it is greater 



2 



than V{xi). As before, the dynamics for the displacements: = Xi — Xi{0), is 
described by the following equations: 



random substrate 



dii/dt = max(0, + 6-1 - 2^) + E- V{i, + i)) 



(2) 



Initially, the particles positions form a ID lattice with unit spacing: XiiQ) = i. 
Pinning forces are local, with a characteristic range A. In our study the pinning 
force is constant inside lattice unit cells, (A = 1), and varies randomly from 
one cell to the next. We chose a binary distribution for V of the type: p{V) — 
p6(y — Vi) + {1 —p)6{V — Vq) , usually with p=0.5. We adopt periodic boundary 
conditions, so that the argument of V in equation (1) is mod{L) + i. 

Both models also depict a discrete one dimensional elastic interface moving 
in the ^ direction in a two dimensional {i,^) plane [10], where the interface 
is initially flat. In the random substrate model a pinning site on the i axis, 
transforms into lines of pinning sites of slope unity in the (i, ^) plane, as 
described in Fig. 2. In the random network case, pinning sites form columns 
of fixed pinning strength Vi in the same plane. 

Approximate scaling properties 

The dependence of the general behavior of the random substrate model on 
the elastic stiffness k, and on the size A of regions where pinning forces do 
not vary, can be simplified by the following scaling relation [19] (provided 
max(«;, kA) < 1): 



Extending this scaling relation to A > 1 results in a hybrid random network- 
random substrate model, so that the pure random substrate model is not 
preserved. However, in the limit A >> 1 the pure random network model 
is approached, so that the random network model and the random substrate 
model approach each other in the very high stiffness limit, as k — > 00. 



3 Dynamical properties 

In the random network model above the threshold Ec all particles eventually 
move with the same finite velocity as a rigid distorted object. Adding the 
N equations of motion Eq. (1) the elastic forces sum-up to zero, so that the 



e(«,A,i) 



K-'^{l,KA,t/K) 



(3) 
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velocity of the center of mass is: [9] 



Urn — E — Ec 



(4) 



where for a chain of N particles Ec — 

The dynamical behavior of the random substrate model is more complex, and 
two distinct regimes appear as the driving force is varied above the depinning 
threshold E^. 

stick — slip regime : Ei < E < E2 ; 

The center of mass of the system moves with a roughly constant velocity, 
modulated by fast fluctuations, although at any time a finite fraction of the 
particles are not in motion. 

The center of mass velocity U obtained numerically for several values of stiff- 
ness K is shown in Fig. 3. Each curve corresponds to a single representative 
random configuration for chains N=100-1000. At the threshold there is a small 
discontinuity, followed by a linear E dependence. The fractions of particles that 
are in motion at a given instant corresponding to the velocity curves in Fig 3. 
appear in the inset. 

Deriving the CM velocity by summing the individual equations Eq. (2) is more 
subtle in the random substrate case, since the equations for particles with zero 
velocity are actually inequalities . Time averages over sufficiently long periods 
are equivalent to those obtained from a linearized version of Eq. (1,2) where 
the ratchet condition is omittted, leading to: 



where n{Vi, E) is the time averaged fraction of particles occupying the i'th 
unit cell. Since the residence time of particles over a strong pinning area is 
longer than over a weak pinning one, the mean occupation fractions niVi, E) 
are monotonously increasing functions on , so that for two systems sharing 
the same distribution of random pinning strengths, the mean velocity of the 
random network is always greater than that of the random substrate system, 
asymptotically approaching it from below in the limit of high velocity. The 
threshold Ei is the solution of the equation: 



URs^E-Y,n{V,,E)V, 



(5) 



E^-Y.n{V,,E^)V, 



(6) 



so that the inequality Ei > Ec is always valid. 
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The center of mass velocity of the random substrate model can be expressed 
as a driving force minus a velocity dependent drag force F{U): 



Analysis of the velocity data shows that F{U) is equal to the static friction 
force El in the limit of zero velocity, increases with increasing velocity above 
El reaching some plateau, eventually tending to E^ for high velocities. This 
weakening of the kinetic friction force with increasing velocity is generated 
by a detachment instability where the fraction of immobile particles abruptly 
vanishes (Fig. 3 inset). Qualitatively similar features were observed in experi- 
ments where an clastic membrane is dragged over a rigid substrate, performing 
stick-slip motion. [7] 

The center of mass velocity for E just above the depinning transition scales as 
[E — EiY . For the random network model, trivially, j3 = 1. Our results for the 
random substrate model also follow a hnear E — Ei dependence, i.e. (3 — 1. 
This is in contrast with Ref. [20], where (3 — 0.47 is reported. The differences 
in exponents may be due to differences between our model and the specific 
realization of the random substrate model in Ref. [20]. This is supported by 
the fact that in nonlinear dynamical systems with quenched randomness, uni- 
versality is often weaker than in equilibrium critical phenomena, so that some 
critical exponents may depend on details on short length scales. A striking 
demonstration of this for a directed polymer in a random medium, a problem 
closely related to the one considered here, was published recently [21]. 

free motion regime: E > E2 ; 

In this regime all particles possess a finite velocity at all times (See Fig. 3, in- 
set). The instantaneous velocity of the center of mass is made up of a constant 
part plus a random washboard like modulation induced by motion over a fluc- 
tuating pinning landscape. In the limit: E » E2, the velocity tends to that 
of the analogous random network Eq. (3), and the relative fluctuating com- 
ponent of the velocity diminishes (Fig. 3). The drag force F{U) approaches 



Our model assumes the presence of local friction forces randomly distributed 
along the interface, while their microscopic origin are outside its scope. The 
stick-slip dynamics in the random substrate model is a consequence of the 
random arrangement of simple local interactions together with the cooperative 
effect of the elastic forces. While the decrease of macroscopic kinetic friction 
with velocity in real materials may have various causes, comparison with our 
results suggests surface inhomogeneity or randomness may play a significant 
role. 




(7) 



Ec. 
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4 Scaling Properties of the Pinned Sate 

For both models, interfaces evolving by the dynamics described in Eq (1) or 
Eq (2), starting from a flat initial configuration, (^j = at t=0) and subject 
to a constant drive E below the pinning threshold, eventually reach a static 
state of strongly strained domains, or strain avalanches [22], seperated by 
virgin particles which never moved. For E approaching E^. from below, the 
random network model undergoes a second order dynamical phase transition 
[9,10]. The following numerical results show that the random substrate model 
follows a qualitatively similar critical behavior, but of a different universality 
class. 

Fig. 2 shows a typical pinned interface for the random substrate model. Here 
interfaces are hindered by the tilted pinning lines, leading to characteristic 
triangular structures with a roughness exponent ( — 1. In contrast, in the 
random network model, distortion of domains is much more pronounced and 
the roughness exponent is = 3/2 [9,10]. Similar to the random network model 
and to many other non-equilibrium systems (e.g. the scaling of avalanche sizes 
in self-organized criticality: [23]), the number of domains of size i per unit 
length in the random substrate model has the generic scaling form: 

n{£) ~ X^'x-^^x) (8) 

where x = £/A(e), and ^(.x) is a slowly varying function for x < 1 with a 
steep cutoff at x > 1. The typical size of the largest domains, A, diverges as 
e~'^ where e = E — E\. The number of domains is not fixed. Close to the 
depinning threshold the number of virgin particles vanishes, and the whole 
system is tiled by strained domains. This global condition implies that the 
first moment of n{l) is a slowly varying function of e that does not become 
singular at criticality, leading to cr = 2 for r > 1, while for r<l,o" = l— r. 
The exponents of the random network model are different: r = 3/2, cr = 2, 
i/ = 2 [9]. 

We obtained n{C,) for the random substrate model by solving Eq (1) numeri- 
cally for an ensemble of 40 random realizations of chains of 10"^ particles. Fig. 
4 shows the rescaled histograms. The collapse of the data for different values 
of e was achieved by rescaling using v = 2 and r = 0.9 ± 0.1. 

The analysis of higher moments of n{i) does not show significant deviations 
from single parameter scaling. The values of critical exponents are siginifi- 
cantly different than those of the random network model, confirming that the 
two models belong to distinct universality classes. 



6 



5 The Phase Diagram 



The properties of the random substrate model depend on the network's elastic 
stiffness. Intuition suggests that when the network becomes very rigid, the 
differences between the two systems should vanish. Fig. 5 represents a {k E) 
phase diagram, obtained from numerical studies of chains of N=100. As the 
stiffness increases, the random substrate threshold Ei approaches E^.. The 
second threshold E2 is of the order of Vmax-, the upper bound of the pinning 
strength distribution, and depends only weakly on n. 

Within the pinned regimes of both models, in region I of the phase diagram, 
configurations of strained domains of both models are roughly the same. The 
dark circles denotes values of E where the relative Hamming distance D [24] 
between locations of strained domains for pairs of systems from each model 
with the same set of pinning strengths {Vi}, is smaller than 0.06. 

For the random network model the largest displacement ^max{^) is bounded 
by A(e)^/^ ~ e^^. The scaling of the random substrate model with stiffness 
given by Eq (3), imphes that for large k > ^^nax the two models track. This is 
a sufficient condition , so it yields a lower bound for the boundary of region I 
where the two models coincide, given by values oi k,E: 

^c-^-^ (9) 



An upper bound for E{k) is given by the relation: E^. — E 1/ k [19]. The 
corresponding boundary is also shown in Fig. 5. The fact that for stiff networks 
the two models are roughly equivalent, means that aspects of friction between 
very stiff solid bodies which are related to inhomogeneities on the interfaces 
can be described by the random network model, or by models based on it, 
elaborated to better conform to real sohds. 



References 



[1] Charge Density Waves in Solids, edited by L. P. Gorkov and G. Gruner (Elsvier, 
Amsterdam, 1989) and references therein. 

[2] H. Fukuyama, and P.A.Lee, Phys.Rev. B 17, 535 (1977); P. A. Lee and T. M. 
Rice, Phys. Rev. B19, 3970 (1970). 

[3] M. Ficgclman, V. Geshkenbein, A. Larkin, and V. Vinokur, Phys. Rev. Lett. 
63, 2303,1989. 

[4] R. Burridge and L. Knopoff, Bull Seis. Soc. Am. 57,341,1967. 



7 



[5] J. M. Carlson and J. S. Langer, Rev. Mod. Phys. 66,658,1994, and references 
therein. 

[6] D. S. Fisher, K. Dahmen, S. Ramanathan and Y. Ben-Zion, Phys. Rev. 
Lett.,78,4885,1997. 

[7] D. P. Villete and J. P. Golub, Phys. Rev. E 47,820,1993. 

[8] D. S. Fisher, Phys. Rev. Lett. 50,1486, (1985); D. S. Fisher, Phys.Rev.B31, 
1396 (1985). 

[9] I. Webman, Phil. Magazine, B56743,(1987); I. Webman, (unpublished preprint, 
1987). 

[10] I.Webman in : Nonlinearities in Complex Systems Ed: S.Puri and S. Dattagupta, 
Narosa, New Delhi, 1997. 

[11] A. A. Middleton and Daniel S. Fisher, Phys. Rev B 47, 3530, (1993); O. Narayan 
and A. A. Middleton, Phys. Rev. B49, 244,1994; O. Narayan and D. S. Fisher, 
Phys. Rev. Lett. 68,3615,1992; Phys. Rev. B46,11520, 1992. 

[12] C. R. Myers and J.P. Sethna, Phys. Rev. B 47,11171,1993. 

[13] G. Parisi, Europhys. Lett. 17, 673 (1992); G.Parisi and L. Pietronero, 
Europhysics Letters, 16, 321 (1991). 

[14] , T. Giamarchi and P. Le Doussal, Phys. Rev. Lett. 72,1530,1994; Phys. Rev. 
Lett. 74, 606,1995. 

[15] O. Pla and F. Nori, Phys. Rev. Lett. 67,919,1991. 

[16] J. Krug, Phys. Rev. Lett. 75,1795,1995. 

[17] L. Balents and M.P.A Fisher, Phys. Rev. Lett, 75,4270,1995. 

[18] In the FLR- Fisher equation is similar to Eq(2), except that there V{i) is 
multiplied by a cosine function of ^. Using the ratchet condition is similar to 
replacing the cosine force by a derivative of a sawtooth potential, except that 
local backword motions down the slope of the sawtooth are excluded. These 
motions do not affect the large scale dynamics since the model has no inertia. 
The critical behavior of the random network model below threshold is in the 
same universality class as the CDW models of Fukuyama Lee and Rice and 
Fisher [2,8], but can be solved analytically. 

[19] L Webman (unpublished) 

[20] D. Cule and T. Hwa, Phys. Rev. Lett. 77, 278,1996; T. Hwa and D. Cule, Phys. 
Rev. Lett. 79, 4930,1997. 

[21] Paolo De Los Rios and Yi-Cheng Zhang, Phys. Rev. Lett. 81,1023,1998. 

[22] The equations of motion Eq(l), Eq(2) have inifinitely many static solution. 
We focus on pinned states reached first from flat initial configurations in the 
presence of time independent external field. 



8 



[23] P. Bak, C. Tang and K. Weisenfeld, Phys. Rev. Lett. 59,381,1987. 

[24] The relevant Hamming distance is: D = ^^^rs (g . ) where is a unit 
step function. 



9 



Random network model 
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Random substrate model 
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Fig. 1. Illustration of the random network and the random substrate models. 
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E/E^ 



Fig. 3. Time averaged center of mass velocity vs. driving force for networks of various 
stiffness. Inset: Mean fraction of particles in motion. 
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Fig. 4. Log-log plot of scaled distribution of sizes of strain avalanches in the pinned 
regime. 
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0.6 1 1.4 1.8 2.2 

Fig. 5. Phase diagram in the (k, E). The network configurations are pinned in phase 
I and phase II. Phase III corresponds to stick-shp dynamics, and phase IV to free 
motion. In phase I the random network and the random substrate models yield 
similar pinned configurations. Theoretical lower and upper bounds of the boundary 
of phase I are shown. 
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